Nonresonant entrainment of detuned oscillators induced by common external noise 
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We have found that a novel type of entrainment occurs in two nonidentical limit cycle oscillators 
subjected to a common external white Gaussian noise. This entrainment is anomalous in the sense 
that the two oscillators have different mean frequencies, where the difference is constant as the noise 
intensity increases, but their phases come to be locked for almost all the time. We present a theory 
and numerical evidence for this phenomenon. 



PACS numbers: 05.45.-a, 05.45.Xt 

Entrainment is a key mechanism for the emergence of 
order and coherence in a variety of physical systems con- 
sisting of oscillatory elements. It is one of the funda- 
mental themes of nonlinear physics to explore the possi- 
ble types of entrainments and clarify their fundamental 
properties. One of the typical entrainment phenomena 
is that caused by an external periodic signal. Consider 
two independent limit cycle oscillators of slightly differ- 
ent natural frequencies u>i, i = 1,2. Suppose that these 
two oscillators are in resonance with the external periodic 
signal, i.e., mtUi — nQ ~ 0, where m and n are integers 
and Q, is the external signal frequency. In this case, it is 
possible that the two oscillators have the same mean fre- 
quency nVL/m. This type of entrainment can be observed 
in various systems as diverse as periodically driven elec- 
trical circuits, lasers with coherent optical injections, and 
biological circadian rhythms. Resonant entrainment can 
be described by using simple dynamical models for phase 
variables and the essential properties are well understood 
(e.g., 0). 

Recent physical and numerical experiments have 
shown that not only a periodic but also a noise-like sig- 
nal can give rise to entrainment between two independent 
oscillators @, H, 0, H, HI • The concept of entrainment of 
limit cycle oscillators induced by common signals has to 
be generalized to include the case of noise-like signals. 
The entrainment by a noise-like signal is a nonresonant 
one in the sense that there is no resonance relation be- 
tween the oscillator and the noise. Entrainment between 
two independent and identical oscillators induced by a 
common noise signal has already been studied 5] : it has 
been analytically shown for a wide class of limit cycle 
oscillators that the phase locking state becomes linearly 
stable by applying an arbitrary weak Gaussian noise. 

However, in real systems, the two oscillators are never 
identical but slightly detuned. We note that the theory 
in Ref. [5J does not guarantee the entrainment between 
nonidentical oscillators at all. As an illustrative exam- 
ple, consider two pairs of oscillators, where one pair has 
a natural frequency lo\ and the other pair has a different 
natural frequency u>2, and suppose that they are sub- 
jected to a weak common noise. The theory tells that 



the phase locking occurs in each pair. However, it does 
not occur between the two pairs because for a weak noise 
the mean frequencies of the pair with lj\ and that with 
W2 are still close to u>\ and u>2, respectively. One might 
expect that the nonidentical oscillators come to have the 
same mean frequency and stable phase locking occurs 
when large enough noise is applied. As we will show, 
this is not the case. It has not yet been clarified at all 
what kind of phenomenon happens between nonidentical 
oscillators. It is necessary to clarify this point for better 
understanding of real systems. 

In this study, we consider a general class of limit cycle 
oscillators and reveal that a novel type of entrainment 
occurs between two nonidentical oscillators subjected to 
a common white Gaussian noise, which we call the non- 
resonant entrainment. This entrainment is anomalous in 
the sense that the two oscillators have different mean fre- 
quencies and the difference is constant even if the noise 
intensity increases but their phases come to be locked for 
almost all the time. 

Let Xi € R N be a state variable vector and consider 
the equation 

X i = F(X i )+6F i (X i ) + G(Xi)ri{t), i = 1,2, (1) 

where F is an unperturbed vector field, 5Fi and 8F2 are 
small deviations from it, G £ R N is a vector function, 
and rj(t) is the white Gaussian noise such that (n(t)) = 
and (r](t)il( s )) = 2DS(t — s), where (■ • ■) denotes aver- 
aging over the realizations of r] and S is Dirac's delta 
function. We call the constant D > the noise inten- 
sity. The noise-free unperturbed system X = F(X) is 
assumed to have a limit cycle with a frequency u). We 
employ the Stratonovich interpretation for the stochas- 
tic differential equation ([I]) . This interpretation allows us 
to apply the phase reduction method to Eq. Jl}, which 
assumes the conventional variable transformations in dif- 
ferential equations. 

If we regard the common noise as a weak perturbation 
to the deterministic oscillators and apply the phase re- 
duction method to Eq. |T]), we obtain the equation for 
the phase variable as follows: 

fa = u + 5uji(<j>i) + Z(&)r;(t), i = l,2, (2) 
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where u> is the frequency of the unperturbed oscillator, 
8u>i is the frequency variation due to 8F il Z is denned 
by Z(fa) = G(X ((j>)) • (gr&d x fa x=XoW ), where is the 
phase variable defined by the unperturbed system X = 
F(X) and Xo(<j>) is its limit cycle solution. By definition, 
Z((f)) is a periodic function, i.e., Z((j)) = Z((j) + 2n). We 
assume that Z is three times continuously differentiable 
and not a constant. It is also assumed that < D/oj <C 1 
to ensure the validity of the phase reduction. 

In order to derive the average equation for fa , we trans- 
late Eq. ([2]) into the equivalent Ito stochastic differential 
equation: 

fa = w + 6uJi(fa) + DZ(fa)Z'(fa) + Z(fa)j](t), (3) 

where the dash denotes differentiation with respect to fa. 
In the Ito equation, unlike in Stratonovich formulation, 
the correlation between fa and rj vanishes. If we subtract 
Eq. © for fa from that for fa and take the ensemble 
average, then we have the average equation 



fa) = (8bJi(fa)) - (Soj 2 (fa)) 

+ D{(Z(fa)Z'(fa)) - (Z(fa)Z'(fa))} , (4) 



where we used the fact (Z(fa)r)(t)) = (Z (fa)) {r)(t)) = 0. 
Each ensemble average on the right hand side can be eval- 
uated by using the steady probability distribution Pi(fa) 
for fa, which can be obtained from the Fokker-Planck 
equation for Eq. ©: i.e., (A(fa)) = A(faPi(fadfa 
where A represents a function of fa. The distribution Pi 
can be obtained as Pi (fa) — l/2n + 0(<Ti, D/u>), where 
Oi = maxo<0<27r \Su)i(fa)/oj\. Since SFi is small, cr, is a 
small parameter. Therefore, Pi can be approximated by 
Pi ~ 1/27T for small D/lu and we have 
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(Z(fa)Z'(fa)) 



2tt 



2tt 



2n 



(Su>i(fa)) ~ — / 5uii(fa)d4> = Suji, (5) 



2n 



Z(faZ'(fad<j> = 0, (6) 



where we used the fact Z(0) = Z(2n). If we substitute 
Eqs. © and © into Eq. ©, we have 



dt 



62) = 5u>i — 5u>2- 



(7) 



Since in general Su>i — Slu2 ^ 0, this equation indi- 
cates that the average phase difference increases or de- 
creases in proportion to the time t. In other words, the 
two oscillators still have different mean frequencies even 
when a common white Gaussian noise is applied, i.e., 
d(4>\) 1 1 dt ^ d(fa)/dt. Intuitively, this result is natural 
because the white noise has a uniform power spectrum 
and does not have a characteristic frequency, which could 
entrain the oscillator frequencies. 

Let 9 and i\> be defined by 6 = fa — fa and i\> = fa + 
fa — 2ojt. The variable 9 measures the phase difference 
between the two oscillators. For small D and Sui.;, it is 



expected that fa still has a mean frequency close to lj. 
Therefore, 9 and ijj can be regarded as slow variables. If 
we change the independent variables form (t,fa,fa) to 
(t, 9, tp) and perform the time-averaging with respect to t, 
we can obtain the Fokker-Planck equation corresponding 
to Eq. © as follows: 

— = -(Su; 1 -6u J2 )—-(6u Jl + Su l2 ) — 

d 



+ D 9^[ 2 {m - r(0)}Q] 
+ D ^[nm + r(0)}Ql 



(8) 



where Q(t,9,?p) is the joint probability distribution and 
r is defined by 



r(6) = ^ JJ z(fa)Z( ( j> + e)dfa 



(9) 



Hereafter we assume the case of 8u>i > Suj2 without loss 
of generality. 

It is in general possible that Z has a period smaller 
than 2-7T. Since Z is not a constant function, we suppose 
that Z(fa) = Z(4> + 2n/n), where n is a positive integer. 
Let h(0) be defined by h(6) = 2{r(0) - r{9)}. It can be 
shown that h(0) > for any 9 6 [0, 2ir). The zero points 
s m of h are given by s m = 2-Km/n, m = 0, 1, . . . ,n — 
1, where so — 0. Equation (jHJ) has the steady solution 
Q s (9) such that it is a continuous function of 9 only and 
satisfies the two conditions (i) Q s (9) = Q s (9 + 2ir) and 
(ii) Jj" Q s (9)d9 — 1. In each interval (s m ,s m +i), the 
solution Q s can be obtained as follows: 



Q s (9) = 



2wh(6) 



exp 



1 



h(y) 



dy 



dx, (10) 



where e — (5lo\ — 80J2) /D > 0. The right hand side of Eq. 
(ITCTjl has singularities at the zero points of h. The value 
of Q s for each s m is given by Q s (s m ) = lim e ^ Sm Q s (9). 
Assume that 9 € (s m , s m +i), i.e., 9 is an arbitrary regular 
point. It can be shown that lim £ _^o Qs(9) = holds 
due to the factor e in the numerator. This implies that 
the probability has to concentrates at the singular points 
s m , m = 0, 1, . . . , n — 1 because Q s satisfies the condition 
(ii). Thus, Q s in the limit e — > is given by 



n— 1 

Q.(9) = -Y l 5(0-s m ), 



(11) 



m— 



where <5 is Dirac's delta function. For small positive s, 
the distribution Q s has narrow and sharp peaks at 9 = 



0,1,. 



1 while Q s is close to zero in the 



regions other than the neighborhoods of these singular 
points. The peaks of Qs become narrower as e approaches 
zero. Equation (fTTjl indicates that multiple peaks exist 
if Z has a period smaller than 2tt, i.e., n > 1. The 
existence of multiple peaks has been pointed out in the 
case of identical oscillators |Q] 
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The above profile of Q s clearly shows that the phase 
locking states, where 9 mod 2ir ~ s m , are achieved for 
a large fraction of time during the time evolution when 
the noise intensity D is relatively large with respect to 
the mean frequency difference 8uj\ — Su>2- i.e., the non- 
resonant entrainment occurs. Let S be a small posi- 
tive constant and Us be the (^-neighborhood defined by 
Us = U^l (s m — S, s m + S), where mod 2tt is taken for 
so — 5. We identify the phase locking state by the condi- 
tion 9 € Us- As shown by Eq. {?]), the present entrain- 
ment is not characterized by coincidence of the mean fre- 
quencies of the two oscillators. Therefore, as a measure 
for the entrainment, we introduce the phase locking time 
ratio fi defined by 



/' 



lim 



Tl 
T ' 



(12) 



where Tl represents the total time length for which 9 £ 
Us happens during the period T. This ratio can also be 
expressed in terms of Q s by /i = Q s (9)d9, where the 
integral is taken over the set Us- Equation (fTTj) shows 
that fj, — ► 1 in the limit e = {5lj\ — 6102) /D — » 0. 

A phase locking state cannot continue for the infinite 
time but phase slips have to happen during the peri- 
ods such that 9 ^ Us because the two mean frequencies 
d(4>i)/dt and d((f>2)/dt are different. Equation j7]) indi- 
cates that the mean frequency difference is given by the 
constant Sloi — 8u>2- This implies that the average num- 
ber of phase slips, which happen in a unit time interval, 
does not become small but remains constant even for rel- 
atively large D compared with 8101 — 5lj2- In other words, 
the average interslip interval remains constant. On the 
other hand, the probability for 9 £ Us decreases and con- 
verges to zero as D increases: i.e., the phases come to be 
locked for almost all the time. These two facts imply that 
a single phase slip completes more rapidly: i.e., the time 
needed for one phase slip decreases and converges to zero 
as D increases. We emphasize that the above mentioned 
behavior is a remarkable feature of the nonresonant en- 
trainment. This behavior is very different from that of 
resonant entrainment by a periodic signal, where the av- 
erage interslip interval diverges and the mean frequencies 
becomes identical as the signal intensity approaches the 
critical value for entrainment. 

In order to demonstrate the above analytical results, 
we show numerical results for an example described by 
the Stratonovich stochastic differential equations 



fa = u>i +siri((f>i)r)(t), £ = 1,2, 



(13) 



where u>i, i = 1,2 are slightly different constants. 

Figure[lja) shows the mean frequency difference Alu = 
d((f>i)/dt — d(cf>2)/dt plotted as a function of D, where io\ 
is fixed to unity and five different values of 0J2 are em- 
ployed. We calculated Alu by replacing the ensemble 
averages in Alu with the time averages, i.e., d(<fii)/dt = 
lhriT^oo{4>i{T) — 4>i(0)}/T. The mean frequency differ- 
ence Aw is not zero except for the case u± = U2 = 1. 
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FIG. 1: Entrainment in phase models with Z = sin(^): (a) 
mean frequency difference Alu vs. D, (b) time evolution of 
phase difference 8 = <j>\ — <j>2 and (c) probability distribution 
Q s (8) for D = 0.02 (x), 0.05 (•), and 0.1 (o), where analytical 
results are shown by solid line. The inset in (c) shows fi 
plotted against D. In (b) and (c), o>i = 1 and 102 = 0.98. 



It is clearly shown that Aw is constant and independent 
of D. This result coincides with the analytical result of 
Eq. ([7]) . The steady distribution Pi((f>i) is approximately 
given by P^fa) ~ (1/2tt)[1 + {D/lwij sin(2<fo) ] for this 
example. This shows that the assumption Pi{4>i) ~ 1/27T 
is reasonable for small D used in the numerical calcula- 
tions. Thus, the result of Eq. (JTJ) holds. 

The time evolution of the phase difference 9 — 4>i ~ <f>2 
is shown for three different values of D in Fig. []Jb), 
where lv\ — 1 and 0J2 — 0.98. These results clearly show 
that the phases are locked near 9 ~ 27rn, ncZ and the 
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FIG. 2: Entrainment in SL oscillators. Probability distribu- 
tion Q s {6) is shown for D = 0.02 (x), 0.05 (•), and 0.1 (o), 
where analytical results are shown by solid line. The inset 
shows Aui vs. D. Parameters are Su>i = and 8u2 = —0.02. 

phase slips occur intermittently. It should be noted that 
the time needed for a single phase slip becomes smaller 
as D increases. This observation is in agreement with 
the analytical result. 

The probability distribution Q s (9) is shown in Fig. 
Hfc) for three different values of D, where u)\ = 1 and 
u>2 = 0.98. The analytical results of Eq. (fTU]) are also 
shown for the corresponding values of e = (u>i — ui2)/D. 
It is seen that Q s is close to the uniform distribution 
for small D or large e. In contrast, the distribution has 
a sharp peak near 8 — for large D or small e. The 
peak in Q s becomes narrower and its position becomes 
closer to 8 = as D increases. This agrees with the pre- 
vious theory since the zero point of h{6) is only 8 = 
in this example and thus the theory tells that Q s has 
a peak only at 8 = 0. It is also seen that the peak is 
not centered at 8 = but shifted to the positive direc- 
tion: i.e., the phase cj>i of the larger natural frequency 
oscillator is kept advanced with respect to <p2 even in the 
phase locking state. The inset of Fig. [ljc) shows that the 
phase locking time ratio fi monotonically increases and 
approaches unity with increasing D. Figure [He) clearly 
demonstrates that the phases are locked for a larger frac- 



tion of the time as D increases. 

In order to validate the theory based on the phase re- 
duction method, we carried out numerical experiments 
for the Stuart-Landau (SL) oscillator 

fa = (1 + icjWj - - V(t), j = 1, 2, (14) 

where ijij 6 C and Cj = 1 + Sujj is a real constant. This 
is reduced to the phase model <pj = 1 + Sujj +sm((j)j)r)(t), 
where <f>j is the appropriately defined phase variable. 

In Fig. [21 the numerically obtained distribution Q s {8) 
is shown for three different values of D, where 8u)\ = 
and 8u>2 = —0.02. The analytical results obtained from 
the corresponding phase model are also shown for the 
corresponding values of e = (Suj\ — 5u>2}/D. A sharp 
peak of Q s appears near 8 = 0. It becomes narrower and 
approaches 8 = as D increases. Agreement between the 
numerical and analytical results is excellent, especially 
in small D region, where the phase reduction method 
gives a good approximation. The inset shows the mean 
frequency difference Acu = d{4>x)/dt — d{4>i)/dt plotted 
as a function of D for the same 5uj\ and 5ui2 ■ It is clearly 
shown that Aw does not depend on D and its constant 
value is given by Alo = Suji — 8u>2- This behavior also 
agrees with the theory. The agreements in the behaviors 
of Q s and Au validate the theory based on the phase 
model. 

In conclusion, we have found the nonresonant entrain- 
ment between two nonidentical limit cycle oscillators 
subjected to a common external white Gaussian noise. 
We theoretically elucidated this phenomenon by using 
a phase model and presented numerical evidence for a 
particular phase model and the SL oscillator. The non- 
resonant entrainment is anomalous in the sense that the 
two oscillators have different mean frequencies, where the 
difference is independent of the noise intensity and the 
average interslip interval is constant, while their phases 
come to be locked for almost all the time for relatively 
large noise. It is expected that a similar entrainment 
occurs for various noise-like signals having broad contin- 
uous power spectra. 
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